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ABSTRACT 

We have modeled the infrared (IR) Spectral Energy Distribution (SED) of cir- 



\Q . cumstellar disks embedded in a HII region and photoevaporated by the external 

ultraviolet radiation. The model applies to the photoevaporated disks (proplyds) 
in the Orion Nebula, most of them illuminated by the 06.5 star ^ 1 Ori C. First 
q ! we calculate the IR emission of a Pre-Main-Sequence star surrounded by a dusty 

globule that is immersed within an HII region. The globule is assumed to be 
^! spherical, homogeneous, optically thin at IR wavelengths and photoevaporated 

according to the Dyson (1968) model. Second, we consider the IR emission of 
a disk directly exposed to the nebular environment. The reprocessing disk is 
passive and treated according to the Chiang and Goldreich (1997, CG97) model. 
We improve over the CG97 treatment by tracing the propagation of the vari- 
ous radiative fluxes (from the star exciting the HII region, nebular, and grazing 
from the disk central star) through the disk superheated atmosphere. Since the 
opposite disk sides receive different amounts of radiation, the flaring angle and 
the surface temperature distributions are different, resulting in well distinguished 
SEDs for the two disk faces. Finally, we combine the globule and disk models 
to estimate the IR emission of proplyds. The energy input from the central star 
and the nebular environment increase the disk flaring angle, and therefore also 
the amount of stellar radiation intercepted by the disk. The relative intensity of 
the disk vs. envelope emission varies with the tilt angle relative to the directions 
of 6 |1 Ori C and of the Earth. We explore the dependence of the SEDs upon the 
tilt angle with respect to the Earth, the distance from 8 1 Oh C, the size on the 
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envelope, the inner disk radius and the temperature of the central star. The 
resulting SEDs are characterized by a broad peak of emission at 30-60 /iin and 
are in general significantly different from those of isolated disks in low-mass star 
forming regions like Taurus- Auriga. Our model indicates that in the presence 
of an external radiation field, relatively evolved "Class 2" objects may display a 
SED peaking at mid-IR and far-IR wavelengths. Also, the model can account 
for the strong mid-IR excess we have recently detected at 10 /im from embedded 
disks in the Orion Nebula. 

Subject headings: infrared: stars — stars: pre-main-sequence — accretion: ac- 
cretion disks — sources: Orion Nebula 



1. Introduction 

Circumstellar disks are thought to be the birthsites of planetary systems, and it is 
primarily for this reason that their properties have been studied extensively since the early 
1980s (Shu, Adams, and Lizano 1987; Beckwith & Sargent 1996; Hartmann 1998). The 
discovery of several tens of extra-solar planetary systems (e.g. Marcy, Cochran & Mayor 
2000) added weight to the argument that disks create planets, although there is some debate 
about how often planetary systems really occur around stars. Nevertheless, an understanding 
of the evolutionary properties of disks around young stars is thought to be an essential part 
of the understanding of our own origins. 

Disks can also affect the early evolution of stars themselves. Several authors have 
suggested that stellar rotation is locked to the orbital rates of the inner disks at early times 
allowing the disk to regulate the stellar angular momentum (Edwards et al. 1993, Choi & 
Herbst 1996), although this idea is still controversial (Stassun et al. 1999, Rebull 2001). 
Typical disk masses and accretion rates appear to be too small to greatly change the mass of 
the stars through accretion (Palla & Stahler, 1999). A number of observable characteristics 
of young star/disk systems play an important role in determining the stellar properties and 
therefore understanding even subtle changes is important to ensure that the stellar properties 
are correctly interpreted. 

Most of the stars in the Galaxy are thought to be born in dense OB clusters such 
as the Orion Nebula (Bally et al. 1998a; McCaughrean & Stauffer 1994). However, our 
understanding of circumstellar disks comes mainly from studies of nearby dark clouds such 
as the Taurus/Auriga complex, and it is not clear that the extant knowledge of circumstellar 
disks applies to the majority of young star/disk systems in the Galaxy. The study of the 
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disks in Orion presents several problems that has slowed progress: they are seen against the 
Orion HII region, a bright, highly nonuniform source of radiation at wavelengths from the 
ultraviolet through millimeter; they appear to be fainter on average at millimeter wavelengths 
(Bally et al. 1998b); and the Orion nebula is three times farther from Earth than the well- 
studied dark clouds. Many of the disks are embedded within the HII region, meaning that 
they are surrounded by ionized halos of gas and dust that complicate the observation of disk 
properties alone. 

To address these difficulties, we present in this paper a series of models of the emission 
from disks embedded within the Orion HII region, under the influence of other radiation 
sources in addition to the stars at their center. Section 2 is a discussion of various sources of 
radiation in the HII region to prepare for the calculations of emission from discrete sources. 
In the following sections, we show how the different radiation affects the IR emission arising 
from three type of sources: a spherically symmetric, dusty globule surrounding a star (section 
3), a star/disk system directly exposed to the ionized environment (section 4), and finally a 
star/disk system surrounded by a dusty globule (section 5). Section 6 contains a discussion 
of our findings with an exploration of the critical parameters, comparison with observations, 
and remarks on the limit of our treatment. 



2. General parameters of the HII region 

We assume the HII region is spherical, isothermal, and has uniform density. For sim- 
plicity, we consider the brightest star at its center as the only source of ionizing photons, 
neglecting the radiation emitted by other stars within or around the nebula. This approx- 
imation is appropriate for the Trapezium cluster, where ^Ori C supplies at least 80% of 
the Lyman continuum luminosity: L^ yc ^8x 10 4 L q versus a total luminosity of ~ 10 5 L Q 
for the entire cluster (Bally et al. 1998a). We assume the HII region is radiation bounded 
with Stromgren radius, Rhii, i-e. all the Lyman-continuum radiation is absorbed within 
the nebula. We use for the HII region the typical parameters of the Orion Nebula, with a 
Stromgren radius Rhii = 1 x 10 18 cm and a distance D = 450 pc. 

To estimate the relative importance of the various dust heating mechanisms, it is appro- 
priate to distinguish between the stellar, F s , and the nebular, F n , radiative fluxes (Panagia 
1974, Natta & Panagia 1976). The stellar flux is comprised of ionizing, Ff uv , and nonion- 
izing, Ff uv , radiation. The former maintains the ionization balance and is mostly absorbed 
in the vicinity of the Stromgren boundary. Once absorbed, it is transformed into nebular 
radiation, F n , that fills the HII region uniformly. F n can be divided into the resonantly 
scattered Ly-a radiation, F^ y ~ a , and the rest, e.g. Balmer continuum and lines, forbidden 
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lines from various atomic species, continuum emission, etc., in short all referred to hereafter 
as nebular radiation, and indicated by F° ther . 

If p denotes the fraction of total luminosity, L S) emitted by the central star in the 
Lyman continuum (Extreme Ultraviolet - EUV; hu > hu LyC = 13.6 eV), the ionizing flux at 
a distance, d, is: 

F s EUV (d)=p^e-^ d \ (1) 

where TEuv(d) is the radial optical depth at the effective wavelength of the ionizing radiation, 
Aw (~ 570A for ^Ori C). The non-ionizing stellar radiation (Far Ultraviolet - FUV; 
6eV< hv < 13.6 eV) is 

Fr v (d) = (l-p)-^ i e-™M (2) 
with an effective wavelength, Xfuv — 1200A (Panagia 1974). 

With the exception of the dusty globules and disks, discussed later in this paper, we 
assume that the amount of dust mixed with the ionized gas is negligible, so that both the 
EUV and FUV fluxes cross the HII region almost unaffected by dust. 

Recombination theory predicts that as long as the electron density is << 10 4 cm~ 3 a 
fraction ~ 2/3 of the Ly-continuum photons are eventually transformed into Ly-a photons 
at \iy~a — 1216 A through the 2p-ls transition, the remaining decaying to Is via two-photon 
emission from the 2s level (e.g. Osterbrock 1989). The generation rate of Ly-a photons is 
therefore 

N Ly-a (3) 

3 hv LyC 

and the corresponding luminosity is 

L L r a = \ P L S ^ \ P L S , (4) 

where \ LyC = c/v LyC = 912 A. 

For the ionizing star we adopt the parameters of an 06.5V star: L s = 6 x 10 38 erg s~ x 
and p = 0.4 (Panagia 1973). 

For the remaining nebular radiation, 

L°: her = pL s - L^ a ^ \ P L S . (5) 

The energy density of the nebular radiation is obtained by multiplying the luminosity 
by the crossing time, Rhii/c, and dividing by the nebular volume, V ne b = {^/^)^R%n- I n 
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the case of the Ly-a radiation, the photons are resonantly scattered many times within the 
nebula before being absorbed. Numerical models (Panagia 1978, Hummer & Kunasz 1980) 
indicate that whereas Ly-a photons undergo approximately 10 5 scatterings before reaching 
the boundary of an HII region, the space crossed between most of the scatterings is very 
short, so that a Ly-a photon reaches the boundary after having traveled on average a path 
length about q ~ 10 — 20 times the radius of the HII region. Therefore, the absorption of 
Ly-a photons by dust is about q times more efficient than it is for non-resonant radiation 
with similar frequency. From the point of view of the dust heating by Ly-a radiation, the 
factor, q, can be treated either as an apparent increase of the Ly-a energy density seen by 
the dust, 

or as an increase of the dust absorption efficiency Q a b s at the Ly-a wavelength, 

Q a6s (Ly-a) = qQ abs (UV). (7) 

This second approach is more correct, since the apparent increase of the energy density holds 
only as long as Ly-a photons are present, i.e. inside a dust-depleted HII region and at the 
outer edges of any embedded globule, disk, etc. In these objects, Ly-a photons will cause 
an increase of the surface temperature, but at an optical depth to UV radiation r ~ 1 they 
will be entirely absorbed. Equation (7) has been adopted in the rest of this paper. 

Concerning the nebular radiation, one has 

Tjother ^PL S ( o\ 



3. Circumstellar globules 
3.1. Neutral flow 

The ionized globules in the central part of the Orion nebula are typically shaped like 
tear-drops with the tail pointing away from ^Ori C. For simplicity, we treat them as spheres 
of neutral gas surrounded by a shell that is photoionized on the hemisphere facing the ionizing 
star. Johnstone, Hollenbach & Bally (1998) and Storzer & Hollenbach (1999) have modeled 
the structure of the globules in terms of photoevaporated winds. They showed that two 
different outflow regimes may develop: EUV-dominated and FUV-dominated flows. The 
former are produced when the ionizing radiation reaches the disk surface, typically in the 
vicinity (d < 10 17 cm) of the ionizing star, and/or in the very early phases of disk exposure to 
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the EUV radiation. FUV-dominated flows are more common and occur when the EUV flux 
is absorbed at large distances from the disk surface. In this case, the disk photoevaporation 
is driven by the FUV radiation. The photodissociation region (PDR) that develops between 
the disk surface and the ionization front may extend hundreds of astronomical units, or a 
few disk radii (~ 100 AU) from the disk. Typical solutions of the flow equations indicate 
that in EUV-dominated flows the density remains approximately constant, whereas in FUV- 
dominated flows the expanding neutral flow has nearly constant velocity until it reaches an 
isothermal shock, which separates the inner region with density oc r~ 2 from the isobaric 
upstream region where the density is approximately constant. 

Storzer & Hollenbach (1999) have shown that in most FUV-dominated cases the shock 
front occurs much closer to the central star than to the ionization front. The region with 
constant density therefore fills most of the globule's volume. The dense inner flow may still 
dominate the radial column density, but it is generated in the outer disk regions, where the 
thermal supersonic neutral wind at 2-6 km s" 1 exceeds the escape velocity from the rotating 
disk. The neutral outflow will be initially cylindrical by symmetry, while a spherical flow 
develops only at distances of the order of the disk size. The star should therefore remain 
in a cavity relatively free of dust, possibly dominated by the presence of collimated jets 
(Bally, O'Dell, & McCaughrean, 2000). In conclusion, both EUV- and FUV- dominated 
flow models predict the presence of broad regions around the central star with uniform 
density. We shall therefore consider only FUV-dominated flows with constant density. Our 
neutral region spans from the ionization radius Ri down the an inner radius i? , encircling 
a dust-free cavity set by the dust evaporation temperature, T e = 1,500K. We shall refer to 
this structure as a "neutral globule" . 

The density within the proplyds is regulated by the dust opacity to the FUV flux (Bally 
et al. 1998a). If the density - and therefore, opacity - increases, the number of FUV photons 
reaching the disk diminishes. The smaller number of FUV photons causes the evaporation 
rate to drop with a consequent decrease of the density (opacity) in the globule. More FUV 
photons will then be able to reach the disk, keeping the system in equilibrium through 
negative feedback. Both Johnstone et al. (1998) and Storzer & Hollenbach (1999) estimated 
the column density of the neutral medium, N(H) ~ 10 21 cm~ 2 , corresponding to Ay ~ 0.5 
assuming average ISM opacity. This column density translates into a radial optical depth 
r Ri ~ 1 — 2 in the UV (both EUV and FUV), the exact value depending on the dust 
properties. In the next section, we shall see how r Ri is directly related to the ionizing flux 
and to the radius of the globule. 
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3.2. Photo-ionized wind 

Behind the ionization front, the gas flow is approximately spherically symmetric. At 
large distances from the globule, the interaction with the wind from the ionizing star will 
create the characteristic stationary arcs of [OIII]+Ho; emission seen around several sources 
in the inner core of the Orion nebula (Bally et al. 1998a, Bally et al. 2000). Dust is trapped 
within the photoevaporating flow, as demonstrated by the arcs closer to ^ 1 Ori C, which 
are prominent at 10/im. In fact, the Ney-Allen nebula is mostly due to the optically thin 
emission of dust in these arcs (Hayward, Houck, & Miles 1994). 

Dyson (1968) first studied the evaporation of dense neutral globules in HII regions 
as a mechanism to provide high turbulence to the ionized medium. In Dyson's model, the 
globules are self gravitating spherical condensations of neutral hydrogen. Under the influence 
of the ionizing radiation, the ionized material at the globule surface streams out at supersonic 
velocities, injecting momentum and energy into the HII region. In this model, the spherically 
symmetric, steady flow has a solution given by 



r / c u 
Ri V v ( r ) 
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rii 

where r is the distance from the center of the globule, n(r) and v (r) are the flow density and 
velocity, and rii is the atomic density at Ri. 

We assume cu = 10 6 cm s" 1 for the sound speed in the ionized gas. The Dyson solution 
neglects the role of gravity. In our case, gravity is negligible at distances from a central 
star r > 100 AU, since the escape velocity, v e = a/ GM±/r ~ 1km s _1 , for a central star of 
mass M+ = 1 M is significantly lower than cu (see also Henney et al. 1996, Johnstone et 
al. 1998). In any case, the emitted IR spectrum is almost insensitive to the details of the 
density and velocity profiles. 

In the expanding ionized flow, the neutral hydrogen and the grains compete for the 
absorption of the ionizing photons. The effects of dust in the photoevaporation flow has 
been studied by Pastor, Canto, & Rodriguez (1991), and we refer the reader to their paper 
for the details. The optical depth within the globule is given by their Equation (12): 

F = n iCn e T ~ + a B n 2 t R t b{T%), (11) 

valid in the plane-parallel approximation, relating the flux of ultraviolet photons entering the 
globule, F, to the opacity from the ionization front to infinity, r^, to r° = a d riiRi. Here, «b 
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is the recombination coefficient to the excited levels of the hydrogen, and a d = 2 x 1CT 21 cm -2 
is the grain cross-section per hydrogen nucleus to the EUV radiation. The function, b(r®), 
has been tabulated by Pastor et al. (1991) for the Dyson's density profile. When r° < 10, 
one can approximate 

6(r°) = 0.12e - 37 ^°. (12) 

In this case (Pastor et al 1991): 

rf = 0.45r d °. (13) 

Figure 1 shows how the optical depth, r°, depends on the ionizing flux (i.e. distance 
from the ionizing star) and on the radius of the ionization front (i.e. globule's size) through 
Equations (11)-(13). In general, r° increases with the globule size and decreases with the 
ionizing flux. At distances from the ionizing star larger than 0.06 — 0.08 pc, it is r° ~ 
1 — 2. This corresponds to the range of values for mentioned above. In fact, under the 
assumption that the neutral, inner part of the globule has constant density, and as long as 
Ro << Ri, r° coincides with the radial optical depth r^. In this case, Equation (11) links 
tr. to the ionizing flux and to the globule radius. Both parameters can be easily obtained, 
or constrained, by the observations. 

3.3. Radiation field inside the globule 

Assuming the number density of globules present within the HII region to be low enough 
that each one of them is directly exposed to the radiation coming from the central star, the 
primary sources of energy within a globule of radius R will be: 

1. The radiation (both ionizing and nonionizing) from the star exciting the HII region: 

W = L ^ (14) 

2. The diffuse radiation, including both Lj-a and other nebular radiation: 

^-Wjsb)' (15) 

3. The radiation from the star at the center of the globule. 

Within the globule, the energy density results from a combination of these three terms, 
plus the thin shell discussed in the previous section, each one attenuated by the dust and 
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gas absorption. Note that now the first term dominates over the second at distances d < 
Rhii/V^P — 0-QRhii- We neglect the IR radiation produced by the globule itself as a 
significant source of dust heating. Figure 2 provides a sketch of the contributions we are 
going to consider, beginning with the ionized wind. We shall assume hereafter R ~ 

3.3.1. Ionized wind 

Figure 3 illustrates how, according to the model described in the previous section, the 
UV flux is absorbed by the gaseous and dusty components within the ionized wind. We 
have assumed d = 0.1 pc, corresponding to ~ 45" at the distance of the Orion Nebula and 
Ri = 100 AU. With this choice of parameters, the FUV radiation is only ps 35% absorbed 
before reaching the ionization front. The FUV radiation also traces the fraction of EUV 
radiation absorbed by the dust 2 . The dust optical depth, depending linearly on the density, 
dominates the absorption of EUV photons for distances larger than r xs 1.3-Rj. When the 
gas optical depth (proportional to n 2 ) becomes dominant, the ionizing flux has already been 
reduced by approximately one half. The fraction of EUV radiation intercepted by the gas is 
therefore predominantly absorbed in a relatively thin layer outside the ionization front. The 
layer acts as a secondary source of recombination radiation, trapped between the neutral 
PDR region inside and the ionized wind outside. Each side receives fifty per cent of the 
recombination radiation, again with Ly-a carrying ~ 1/3 of the energy and FUV and optical 
lines carrying the rest. A crude estimate indicates that the flux emitted in each direction is 
approximately 1/2 x 1/2 x 1/2 = 1/8 of the flux directly coming from the ionizing star, the 
three factors accounting for the fraction of the energy emitted in the EUV, the fraction of 
EUV radiation absorbed by the shell, and the fraction of photons emitted up/downstream. 
This amount is not negligible and has been included in the model. 

We treat the radiative flux within the ionized wind in a 1-d approximation. The flux, 
F s , coming from the ionizing star has a cylindrical symmetry along the globule axis, defined 
by the line joining the globule center to the ionizing star and equations (1) and (2) can 
be used with the optical depth, r^°, introduced in Section 3.2. The nebular radiation flux, 
F n , as well as the radiation produced at the ionization front, Fif, are regarded as uniform, 
extended sources at the outer and inner edge of the ionized region, respectively. Treating 
the dusty region as a plane parallel slab illuminated by two uniform radiation fields, we use 



2 In this paper we assume Q v = 1 for the dust absorption efficiency in the EUV, FUV, Ly-a, and nebular 
radiation. The corresponding optical depths are therefore equal. 
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the attenuation averaged over 2tt steradian: 

1 

<e^> 27r = J e-^dx (16) 
o 

In general, all nebular Lj-a photons entering the ionized wind will be rapidly absorbed by 
the dust, whereas a fraction, (e~ r ^ c ) 2 , of the nebular radiation will cross the region. 

To these terms, we add, together with the flux from the central star, F*, the fraction 
of nebular recombination lines, Ly-a excluded, coming from the back side after crossing the 
neutral globule. 



3.3.2. Neutral globule 



Within the neutral region, we use a 2-d treatement. No EUV radiation is present, and 
the fraction of FUV radiation crossing the ionization front is further attenuated by an optical 
depth, 

2 



T 2 d(r, 9) = T Ri 



14 (£) -M^). B h ) 



1/2 



;i7) 



between the point Q{9 h ) at the globule surface and the point P(r, 9), along the line parallel 
to the globule axis. It is 9 h = arcsin [r sin 9 / R t ] (Figure 4). Note that since the globule is 
homogeneous, distances and optical depths are related by the simple relation: 



Ri 



(18) 



On the hemisphere facing the ionizing star, the ionized shell converts UV photons into 
recombination line photons. Like the ionized wind, the neutral globule is optically thick to 
Ly-a photons and translucent to recombination lines at the range of optical depths we are 
dealing with. To illustrate their propagation, we consider first the energy input, E, from an 
isotropic radiation field with specific intensity, I u , on the surface of radius, Rf. 



E v = (A7iRf)nI u . 
Within the globule, the radiation field is attenuated by a factor: 



-T'(r) 



sphere 



1 [* P -T'(r,6>) 

2 Jo 



sin 9'd9\ 



(19) 



(20) 



where r'(r, 9') is the optical depth between the point R, assumed for simplicity on the globule 
axis (9 = 0) and a point S(9') on the surface (see Figure 4). Figure 5 shows that when the 
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globule is optically thick, the energy density at the center goes to zero, whereas at the surface 
it approximates 50% of its value in the open space, as the flux comes only from the outer 
side. Continuity across the globule edge is preserved, as the outer energy density is also 
reduced in the globule shadow with respect to the value in open space. The attenuation in 
the vicinity of the surface can then be estimated integrating Equation (20) between and 
tt/2 in the limit >> r, r being the radial optical depth from the surface. In this case, 
Equation (20) reduces to Equation (16). 

Since the two hemispheres are illuminated by different sources of radiation, isotropy is 
broken, and Equation (20) must be replaced by an equation weighting the two contributions 
separately. The flux at a point P will now depend also on the angle 9 (see Figure 4). It is 
convenient to write Equation (19) as 

na 2 (4nl u ), (21) 



and split the isotropic flux intercepted by the sphere cross section in two terms: 

4nl„ = ilMl + n 2 A 2 I 2 u . (22) 

Here, Qi are the solid angles subtended by the areas radiating with intensity, Ii, and Ai are 
the corresponding mean attenuations. It is 

»27r r 0'^{rfi,4>') 



and, 
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A 2 (r,9) 
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(23) 



(24) 



The optical depth, r'(r, 9, 9', <fi'), between the point, P(r,9), inside the globule and a point, 
S(9', 0'), on the surface of the sphere is related through Equation (18) to r', solution of the 
equation 

r' 2 + 2r (cos 9 cos 9' + sin 9 sin 9' cos 0') r' + r 2 - R 2 = 0, (25) 

whereas the angle, 9' 27T {r,9,(f)'), at P between the direction of the ionizing star, and the 
globule's equator is given by the equation: 

2r 2 sin 9 cos 9 cos d> tan 9L + r 2 sin 2 9 - R 2 = 0. 



r 2 cos 2 9 tan 2 9' 2n 



(26) 



3. 3. 3. Flux from the central star 



The flux from the central star depends on r 2 and on the radial extinction. Since 
the central stars of the proplyds in the Orion Nebula typically have late spectral types, 
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we account for the corresponding reduction of the dust optical depth through the ratio of 
absorption efficiencies Q^XtJ/Quv — Q(^t+), where A^ is the effective wavelength of the 
photospheric radiation. At distances larger than R , the stellar flux is 

Fjr) = J^-e~ Trir)QiXT *\ (27) 

where r r is the radial optical depth from the central star. Figure 6 shows the flux densities 
along the symmetry axis for a globule again with i?j = 100 AU and d = 0.1 pc. For the star, 
we assumed T* = 4,000K and R* = 2.5i? . We shall refer hereafter to these values as the 
standard case. Within the ionized wind, the dust heating is dominated by the FUV flux. 
Inside the globule, the dominant source of heating is the radiation from the central star, in 
the inner parts, and the Lj-a radiation closer to the ionization front. 



3.4. Infrared emission 

The dust temperature, Td(r, 9), is determined by the local energy balance of dust grains. 
Within the neutral part of the globule, it is given by the equation 

POO 

o 9 {Quv [F FUV (r, 9) + F Lya (r, 9) + F REC (r, 9)] + Q(A T jF*(r)} = Aa g / Q(u)nB [u, T d (r, 9)] du, 

Jo 

(28) 

where Frec refers to the recombination radiation emitted at the ionization front and a g = 
na g 2 indicates the geometrical cross section of a typical grain of radius a g . A similar equation 
holding for the ionized flow includes the EUV flux, and, in our approximation, has the 
distance, r, as the only independent variable. 

The dust temperature depends strongly on the assumed grain parameters. Following 
Natta & Panagia (1976), we may take for reference "type-1" grains with a g = a\ = 0.1 //m 
and emissivity efficiency in the infrared 

QM = ^, (29) 

and the smaller "type-2" grains with a g = ai = 0.02 /iin and emissivity efficiency 

Q[y) = 7 x 10~ 4 ^|. (30) 

Type-1 grains have properties similar to the classical Draine and Lee (1984) silicates. Power- 
law emissivities allow analytical solutions to Equation (28) for T d . Writing the last two 
equations as 

Q{y) = Q A _/3 , (31) 
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the black body averaged dust opacity Qp is 



Q^= / Q(u)7rB(u,T d )du 



1.510 x 10- 4 Q TI 
5.150 x HT 4 QoT d 6 



= 2. 



1; 



(32) 



Equation (28) therefore becomes 



[F FUV (r, 9) + F Lya (r, 9) + F REC (r, 9)} + Q(A r jF*(r) 

4-Q/3 



V(4+/3) 



T d (r,0) = 



(33) 



and a similar equation holds for the ionized wind. In Figure 7, we compare the temperature 
distributions along the globule axis resulting from the two types of grains, assuming the 
standard parameters. The smaller grains reach temperatures more than twice that of the 
larger grains, i.e ~ 260 K vs. ~ 120 K. Also, with the exception of the region closer to 
the central star (r < l/3i?j), the temperature is almost constant due to the weak dependency 
of the temperature on the flux. In the rest of this paper, we shall use Type-1 grains to be 
consistent with our treatment of the disk, as also assumed by Chiang and Goldreich (1997, 
CG97). This choice allows us to use Equation (33) to obtain the dust temperature and to 
apply the tabulated dust efficiency of Draine and Lee (1984) to calculate a more realistic 
spectral energy distribution. We remark here that scattering of stellar radiation from dust 
grains is not considered in our model. 

The emitted IR spectrum is given by the integrals: 



for the ionized wind. To keep our treatment as simple as possible, we neglect the absorption 
of IR radiation in the energy budget of the envelope. This assumption is valid within the 
range of the optical depths in the UV and allows us to use Equation (33) avoiding the iterative 
computation of the IR radiative transfer. On the other hand, the 10 /iin silicate feature can 
be self-absorbed even at our relatively modest optical depths, having an absorption cross 
section comparable to that at 2.2 jum. Therefore, we introduced the exponential factors in 
Equations (34)- (35) to account for this effect. 

Figure 8 shows the Spectral Energy Distributions (SEDs), vF IR (is), calculated for the 
Draine and Lee (1984) dust grains and, for comparison, for the Panagia's type 2 grains. 
The resulting SEDs peak at ps 25/zmand 10/xm, respectively. Depending on the grain 
properties, the emission can be dominated by neutral globule or by the ionized shell. 
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in the neutral part of the globule, and 




(35) 



4. Externally illuminated disks 



Our next step is to evaluate the IR SED of circumstellar disks embedded in HII regions. 
Although standard models for T-Tauri disks assume the star+disk systems to be isolated, 
a number of studies have been devoted to the effects of the surrounding environment on 
the disk thermal structure and emitted spectrum. In particular, Natta (1993) and Butner, 
Natta, & Evans (1994) modeled the disks' thermal emission when they are surrounded by 
envelopes. Both studies find that at large distances from the central star the disks become 
significantly warmer than if direct heating alone was considered. The heating mostly affects 
the long wavelength emission from the disk. One may expect to obtain similar results for 
the disks in an HII region. 

As for the globules, the disk temperature depends on a combination of heating by the 
central star in the disk, external heating due to the radiation directly coming from the 
ionizing star, and heating from the nebular environment. We neglect heating from accretion, 
which should be a minor component outside a few stellar radii. 

The radiation from the nebular environment typically hits the disk surface at angles 
much larger than a, the grazing angle of the radiation coming from the central star. In 
particular, the flux coming directly from the ionizing star depends on the tilt angle, d , of 
the disk axis with respect to the direction of the ionizing source and acts only on one of 
the two disk faces. The nebular radiation, isotropic and roughly constant within the nebula, 
flows uniformly into the entire disk surface. For a geometrically thin, optically thick disk it 
is a 0. AR^/a (Ruden & Pollack 1991) and the energy input from the disk's central star 
will be dominant over the external radiation at disk radii 
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on the opposite face, having used the standard parameters and having indicated with 
d$.i = <i(pc)/0.1 pc. In Figure 9, a d is plotted as a function of the disk tilt angle d for 
various distances from the ionizing star. The dashed line represents the limit of pure nebular 
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illumination. This occurs a) on the disk face opposite to the ionizing star, b) on both faces if 
the disk is oriented edge-on with respect to the ionizing star, (9 d ~ 90°), 3) at large distances 
from the ionizing star. Figure 9 clearly shows that, with our standard set of parameters, the 
heating of typical circumstellar disks is generally dominated by the nebular environment at 
distances larger than a few AU. 

In reality, the EUV flux hardly reaches the disk surface in environments like the core of 
the Orion Nebula, since the neutral photoevaporated wind pushes away the ionization front 
from the disk. All radiation fields, including the radiation from the disk's central star, are 
attenuated by the evaporated wind. In the next section, devoted to the actual predictions for 
the SEDs of photoionized globules in the Orion Nebula, this effect will be explicitely taken 
into account. Here, however, we shall not reduce the external fluxes by generic absorption 
coefficients. This because, on the one hand, the main goal of this section is to illustrate how 
we build the model, and we want to keep the formalism as close as possible to our assumed 
reference model for circumstellar disks. On the other hand, at large distances from the 
ionizing star, photoevaporation can be so negligible that the diluted EUV flux can actually 
reach the disk surface. Prominent photoevaporated proplyds in the Orion Nebula are clearly 
concentrated within pa 30" from ^ 1 Ori C. In particular, one does not expect the rear face 
of the disks to be photoevaporated, as it is heated only by the disk's central star, and by 
the nebular radiation possibly attenuated by the evaporated material streaming away from 
the system (the cometary tails of the proplyd). The dark disks seen in silhouette against 
the nebular background do not show evidence for photoevaporeted winds (McCaughrean & 
Stauffer 1994). In principle, they can can be modeled using the equations presented in this 
section. The reader is cautioned, in any case, that a comparison with real sources should 
take into account absorption effects 3 . 



4.1. Disk Spectral Energy Distributions 

We calculate the disk emission using our modified version of the CG97 model, which 
provides a self-consistent treatment of the vertical hydrostatic equilibrium and radiative 
transfer of a disk. Disk flaring results from the hydrostatic equilibrium of the disk material 
in the gravitational field of the central star (Shakura & Sunyaev 1973, Kenyon & Hartmann 
1987). In the CG97 model, radiation from the central star grazes the disk surface at an 
angle a and is entirely absorbed once it has reached a tangential optical depth tii ~ 1 at 



3 For reference, the envelope discussed in the previous section reduces the flux from the ionizing star by 
a factor w 1/5, an amount comparable to moving the source twice as far. 
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the characteristic wavelength of the stellar radiation. This absorption produces an optically 
thin layer of "superheated" dust grains, since the grain emissivity in the infrared is lower 
than the absorption efficiency in the visible. Half of the radiation emitted by the layer goes 
inward and regulates the inner disk temperature. 

We have modified this scheme to account for the external radiation field. Since the 
external radiation usually hits the disk at an angle much larger than a, photons entering 
along the normal to the disk surface will be absorbed deeper in the disk than those coming 
from the central star. We introduce a perpendicular optical depth r_|_ related to the tangential 
optical depth by the equation r± t\\ol. Instead of referring to the ry = 1 limit for the grazing 
radiation as CG97, we use the optical depth, r± as the variable that governs the propagation 
of the radiation through the disk. In this way, the dust temperature can be estimated at any 
disk radius, a, as a function of rj_ instead of the single value estimated at r = 1 by CG97. 

With these assumptions, the flux from the central star is attenuated by a factor, 
e -T ± Q(A T j/a(a)^ w ithin the disk, whereas the flux from the ionizing star is reduced by a 
factor, e - T ±/ cos9 d^ jo be consistent with the notation used in the previous sections, t_i_ is 
now the optical depth at UV wavelengths. In principle, the disk flaring produces a mod- 
ulation of 9d, but we ignore this effect considering that the variations average to zero over 
the circumference. We also assume that all the radiation coming from the ionizing star is 
locally absorbed by the dust, directly or through absorption of recombination radiation. It 
is, therefore, 

F^) = ^ 2 cos(6 d ).e-^^. (40) 

The attenuation of the nebular radiation within the disk can be readily expressed as a 
function of r± by using Equation (16), assuming the disk curvature is negligible: 
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where x is the distance from the normal to the disk surface. The two terms on the right side 
represent the Ly-a and the nebular lines, respectively. The total external flux on the disk 
face pointing to the ionizing star is therefore 

F ut(r ± ) = F*(t ± ) + F nebular (r ± ), (42) 

whereas on the other face only the term containing F ne f m i ar (T±) will be present. 

The external radiation affects the geometry of the disk. In particular, the grazing angle, 
a(a), depends on the external flux. Following CG97, it is: 

„(„) (43) 
a da \a J 
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where a,\ ~ OAR+, and H is the disk height above the midplane. For a disk vertically 
isothermal, it is 

?"(*)"£)"- 

where the effective temperature, T e , is a measure of the energy input at the disk surface and 
T c = GM*\xj (kR+), n being the mean molecular weight of the gas. If an outer radiation field 
is present, it is 
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With this expression for T e , Equation (43) cannot be solved analytically. In the limit of outer 
radiation field dominating the disk heating (i.e. at large distances from the disk central star), 
the first term can be neglected and we have: 
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In this case, 
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In Figure 10, we show the vertical disk profiles calculated for a disk facing the ionizing star 
at d — 0.1 pc, with no extinction in front of the disk surface. In this situation, the flaring 
calculated on both faces, i.e. with and without the EUV radiation, is significantly higher 
than that resulting from the standard CG97 model. In practice, these disk profiles are only 
indicative of an optical depth surface, since when photoevaporation is present there is no 
discontinuity between the disk surface and the envelope. A large flaring angle can rather be 
regarded as an indication for disk photoevaporation. The case with no direct radiation from 
the ionizing star provides the minimum flaring profile for the back surface of a disk within 
the HII region (assuming negligible extinction from the envelope material). 



The disk flaring will be dominated by the external flux at distances larger than 
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With our typical parameters, R flaring — 2AU. 
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The temperature of the superheated dust layer is given by the combination of internal 
and external contributions at the optical depth, t±: 



T ds (r±,a) 
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is a characteristic temperature of the dust grains, related to the black body averaged dust 
emissivity at temperature T by 
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In Figure 11, we plot the temperature distribution calculated at various distances from 
the central star as a function of the vertical optical depth. By comparison, the CG97 
approximation of uniform temperature across the superheated layer: 
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(53) 



would have provided T ds = 355 K, 220 K, 147 K, and 113 K for the a = 3, 10, 30, 100 AU, 
respectively. 



The flux emitted by the superheated layer will be given by 
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having assumed that all the flux is absorbed when t± = 3. The luminosity density from the 
superheated layer, normalized to the stellar luminosity L+, is 
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where Oj = 0.07 AU and a Q = 100 AU are our assumed inner and outer disk radii, and we 
have taken into account that only 1/2 of the flux is radiated in the outer space. 

In what concerns the inner disk, we adopt the same scheme of CG97 with modifications 
needed to account for the ambient radiation. Equations 12a-c of GC97 now become: 
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having defined C4 as 

c 4 =(f) v Vr.. (57) 

The three regions are limited by the radii a a \ b and a b \ c , derived by setting r_|_ = 2 and 
r± = 1.5, respectively. With CG97, we have defined 

r(a) = Qi, 2 /€yS a- 3/2 (58) 

with Qi = (^jF^j and Q 2 = (jf^j f° r the two conditions, Ky = 400 cm 2 g -1 , and S = 
10 3 gcm -2 . Consistent with CG97, we used Equation (53) for the temperature, T^ s . With 
our standard parameters, a a \ b = 140 AU. The temperature of the internal disk is constant at 
Tdi — 60 K over most of the disk. The luminosity density from the disk interior normalized 
to the stellar luminosity L± is now 

L in (y) = ^ r [l - e~ T ^] B [u, T in (x)\ xdx. (59) 
L * Jai 



Figure 12 shows the SED for the front surface of the flared disk. The distribution reaches 
a minimum at ~ 4 /im, and then rises up to a maximum at ~ 60 fim. The silicate peaks 
prominent at 10 and 20 fj,m, are emitted in the superheated layer. In Figure 13, we present 
the SED for the face of the disk opposite to the ionizing star. There is a general reduction of 
the disk luminosity, both from the superheated atmosphere at short wavelengths and from 
the interior in the far-IR. As we already mentioned, this type of SED can be representative 
of the dark silhouettes seen in projection against the nebular background in Orion. For 
these disks, Ly-a radiation may play a role only if the disks are still embedded within 
the HII region, whereas the other nebular radiation is relevant also if the disk is in the 
PDR region behind the ionization front. One must remember that the Orion nebula is 
essentially density bounded in the direction of the Earth, being limited by a nonuniform 
veil of foreground material which is not completely optically thick (O'Dell & Yusuf-Zadeh, 
2000). Dark silhouette disks can therefore be subject to a broad range of radiative input, as 
shown for example by the presence of [OI] emission in only one object reported by Bally et 
al. (2000). 



5. Photo-ionized proplyds 

In this section, we model the thermal emission produced by a photoevaporated disk 
enshrouded in a dusty photoevaporated envelope. We treat the envelope as a dusty pho- 
toionized globule of the type discussed in Section 3, whereas the disk treatment will follow 
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the model presented in the previous section. However, when the globule is merged with the 
disk, some of the assumptions adopted for the isolated cases are no longer valid: 

1. The disk is now directly exposed to only a fraction of the outer radiation, due to the 
globule attenuation. In particular, no EUV radiation reaches the surface of the disk. 

2. The disk is exposed to the radiation reemitted by the globule, mostly in the infrared. 
We shall assume that the IR radiation crosses undisturbed the superheated layer and 
is absorbed in the inner disk. 

3. The feedback contribution of the disk emission to the heating of the globule can be 
neglected, since a passive disk emits mostly in the far-IR, where the dust opacity is 
negligible. 

4. The radiation emitted by the disk and by the central star is attenuated crossing the 
envelope. We also take into account the attenuation of the envelope radiation within 
the envelope itself. 

We shall assume the globule to be spherical with the disk lying in the equatorial plane. 
The angles, 0d, between the disk axis and the direction of the ionizing star, and, # e , between 
the disk axis and the Earth are both critical to model SEDs matching the observational 
data. They have been taken into account in an approximate way: the former modulates 
the amount of stellar radiation falling on the disk surface through Equation (40) but not on 
the envelope. The latter affects the tilt angle of the disk with respect to the observer. For 
simplicity, however, the plots presented in this section have been obtained assuming both 
9 d = 0° and 6> ffi = 0°. A disk oriented face-on both to the ionizing star and to the Earth is 
located along the line of sight of, and behind, the ionizing star. 

The globule structure is described by the same set of equations discussed in Section (3). 
The only change is that the term containing 1% in Equation (22) is now equal to zero, 
since the disk, optically thick to ultraviolet radiation, shields the hemisphere opposite to the 
ionizing star. In what concerns the radiation flowing into the disk, both the fluxes coming 
from the HII region and from the central star are attenuated by the globule. In agreement 
with our previous assumption of FUV dominated flows, the EUV flux is entirely absorbed 
before reaching the surface of the disk. The radiation propagates through the superheated 
layer and then through the disk interior. Only the IR flux re-emitted by the globule crosses 
the thin layer undisturbed. 

The total IR flux emitted by the globule and reaching the disk surface is given by 

F(a)= / / / / em, V ' cos 9 a r sin OdudrdOd^ (60) 

Jo Jo Jr % Jo 4 ™a 
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where F env = nd^a 2 g {y)B{y-,T l i) is the flux emitted by a cubic centimeter of dust within the 
globule (cfr. Equations (34-35)), and d a is the distance between the point, P(r, 0, 0), within 
the envelope and a point, A, on the surface of the disk at a projected (equatorial) distance, 
a, form the star and on the plane defined by = 0. It is 

d a = \Jr 2 + (2a sin cos — cos 9 tan a) r + a 2 (1 + tana)), (61) 

and a is the angle between the line, PA, and the normal to the disk surface, given by: 

r sin 9 cos r cos — a tan a 

a = sin a : h cos a (62) 

d a d a 

An accurate estimate of d a would require an iterative solution of Equation (60), since the 
actual position of the point, A, within the globule depends on the flaring angle, a (a), which 
in turn depends on the radiative heating on the disk. To avoid the rather heavy numer- 
ical computation, we adopted the approximate assumption that the position of point, A, 
is initially given by Equation (47), this time with the absorption through the envelope ex- 
plicetely taken into account, and then used the value of F(a) resulting from Equation (60) 
to recalculate the disk flaring angle. 

The flux emitted by the disk is calculated following the treatment presented in Section 4 
with a few changes. We assume that the temperature of the superheated layer does not 
depend on the flux emitted by the globule, as the IR emission propagates almost undisturbed 
through the disk atmosphere. Due to presence of the absorbing envelope, the flaring angle 
now depends on the distance, a, also via the external flux, F ou t, (Equation 42). In general, 
Rfiaring (Equation 49) remains of the order of a few AUs over a broad range of parameters, 
and even when the extinction by the wind/envelope is included, as it affects both the flux 
from the disk central star and from the outer space. Within the Orion Nebula the disk 
flaring is usually dominated by the external flux, so we shall adopt the approximation of 
Equation (47). The luminosity of the superheated layer is given by Equations (54)-(55). 
In what concerns the temperature of the disk interior, it will depend on the sum of the IR 
radiation received from the superheated layer and from the globule, plus the flux emitted 
by the disk interior itself when the disk is optically thick. Equations (56)- (59) apply also in 
this case. 

Neglecting the interstellar contribution, the optical depth between P and an observer 
on Earth is 



r„(r, 0, 0, 0e) = T R,Q{ V ) 
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(63) 



The angles, 9 R , and, 0/?, define the point where the line from P with angle, 0®, intercepts 
the globule surface. They are related to the coordinates of P and to e by equations: 

r sin sin = i?; sin 0# sin(7r — 4> R ) (64) 
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and 

rsin0cos0 + -R i sin0 Ji cos(7r-0 Ji ) 

tan# ffi = 5 7 7 65 

itj cos or — r COS 

solved numerically. Crossing the ionized wind, the radiation is further attenuated by an 
optical depth ~ Qv^T • We neglect the self absorption of the IR radiation emitted by the 
ionized wind 

Figure 14 shows the SED distribution arising from the front face of the disk. With our 
standard set of parameters, the emission is dominated by the disk atmosphere at wavelengths 
< 40 fim and by the internal disk emission at longer wavelengths. The envelope emission 
never produces a noticeable contribution. The disk atmosphere has a pronounced peak of 
emission at 3/zm due to secondary heating by the envelope. The relevance of this effect 
depends on the amount of dust located in the immediate surroundings of the star. We have 
assumed uniform density within the globule up to the dust evaporation radius, which is 
of the order of a few stellar radii. A larger dust free cavity would significantly reduce the 
intensity of the 3 /zm bump. 



6. Discussion 



The geometry we have assumed to illustrate the physical basis of our model represents 
in practice an extreme case. A disk oriented face-on with respect to the ionizing star receives 
the maximum amount of radiation. This causes a substantial increase of the flaring angle, 
and therefore of the inner stellar flux intercepted by the disk. If the disk is also oriented 
face-on with respect to the Earth, the observer will receive the highest amount of disk 
radiation. Even if chances of having such a fortunate alignment are rather low, one of the 
most intriguing 10 /im sources of the Orion Nebula, SC3 (McCaughrean & Gezari 1991, 
Hayward et al. 1994; Figure 15) can be readily interpreted in this way. With a 10 fxm flux 
of ~ 4Jy, SC3 is the brightest compact source in the Orion Nebula, excepting the BN/KL 
region. It is located at a projected distance ~ 1"W of ^ 1 Ori C, and mostly for this reason 
it has never been observed at wavelengths shorter than 2 /im. In our 10 /zm images, SC3 is 
spatially resolved with FWHM ~ 0.5", corresponding to approximately 200 AU, the typical 
size of circumstellar disks in Orion. A direct comparison with the nearby complex of arc-like 
structures produced by the wind-wind interactions suggests that SC3 is a different kind of 
phenomenon. Despite the fact that it has the shortest projected distance from ^Ori C, 
SC3 is apparently undisturbed by the powerful wind from the 06.5 star, as it is the only 
extended source with circular shape. This indicates that the physical distance between 
SC3 and ^ 1 Ori C is much larger than for the other sources. The simplest interpretation 
is that SC3 is a star+disk system in the background with respect to 9 1 Ori C with a very 
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high mid-IR brightness produced by the disk being oriented face-on both to ^Ori C and 
to us. Figure 18 shows that the 10 /im photometry of SC3 is compatible with our standard 
model and a distance slightly lower than 30" =0.065 pc. It must also be noticed that the 
low resolution 10 /im spectrum of SC3 obtained by Hayward et al. (1994) shows silicate 
emission, whereas our standard model predicts at these distances an almost featurless SED 
dominated by the inner disk emission. Further data on this sources are needed to constrain 
our set of parameters. 

In general, disks are randomly tilted with respect both to the ionizing star and to the 
Earth, and therefore their apparent brightness will be reduced in comparison to our standard 
case. In Section 4.1 we stated that the tilt angle with respect to the ionizing star, 0^, can be 
simply accounted for by a projection factor cosOj, in the flux received by the disk surface. 
However, this is strictly true only for tilt angles 6a < a(a ). Since disks are intrinsically 
flared, large tilt angles may cause part of the surface exposed to the ionizing star to go under 
the edge's shadow. When this happens, the heating of the surface is no longer uniform, the 
flaring angle varies with the azimuthal angle and the disk loses its symmetry around the 
rotation axis. Along its orbit around the central star, the disk material will experience a 
periodic modulation of pressure and temperature, and the resulting differential drag may 
cause the disk to warp. Photoevaporation will occur only from the area directly exposed to 
the flux from the ionizing star, as well as on the opposite disk edge. Photoevaporation from 
the edge appears compatible with the observations showing that disks in the Orion Nebula 
are truncated (McCaughrean, Stapelfeldt & Close 2000). As the analysis of the structure 
and evolution of circumstellar disks in these conditions goes beyond the scope of this paper, 
one can maintain the "cos^" treatment to account for a reduced flux falling on the disk 
surface, keeping in mind that the derived tilt angle in this case is an overestimate of the real 
one. 

In what concerns the tilt angle with respect to the Earth, e , Chiang and Goldreich 
(1999) have studied the corresponding variation of the SEDs for isolated disks. They found 
that the SED shows negligible variations until 6® approaches the flaring angle at the disk 
outer radius ct(ao). When the inner disk regions enter in the shadow, the short wavelength 
emission begins to disappear from the vantage point of the observer. To account for this 
effect, we simply multiply the flux emitted by the disk at each radius by the fraction of visible 
area. Figure 16 shows the dramatic reduction of the disk brightness at short wavelengths, 
together with the stellar radiation, when 6> ffi increases from 45° to 50° crossing a(a ) ~ 48°. 
Our geometrical treatment is only valid as long as the disk is optically thick at all wavelengths 
and large enough to contain a strongly flared surface. In the scheme of CG97, the outer disk 
regions are in a completely optically thin regime (region c of Equation 56), and therefore 
the transition is modulated by the variation of optical depth across the disk. In our case, 
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however, the truncated disk remains optically thick (region a) up to the disk edges. We find 
that the disk remains optically thick up to a Q ~ 100 AU over a wide range of parameters. In 
particular, our standard disk+envelope model has a a \ b = 113 AU. Thus, the change of SED 
at short wavelengths is much more abrupt. The decrease of disk emission at A > 30 /im, is 
correct as long as the outer parts of the disk are also optically thick and, therefore, sensitive 
to the disk orientation. In conclusion, the disk orientation with respect to the Earth can be 
regarded like a "switch" for the short wavelength radiation. If the star is directly visible (e.g. 
with HST), then the brightness at near- and mid-IR wavelengths will be high, otherwise if the 
star is invisible. Transition objects with heavily reddened stars should be rare. Observations 
at 10 fim of dark silhouette disks appear to confirm this result: the the large dark silhouette 
114-426, seen edge-on, has a star under A v > 60 magnitudes of extinction (McCaughrean et 
al 1998) and remained undetected at [N] ~ 9 in our deep survey performed with the UKIRT 
3.8 m telescope (Robberto et al, 2002), whereas the dark silhouette with the brightest central 
star, object 218-354, has been detected at 10 /im both by Hayward & McCaughrean (1997) 
and by us. 

In the standard case, the envelope contribution is negligible at all wavelengths with 
respect to that coming from the disk. For some disk orientations, however, the envelope 
emission becomes important. This because the envelopes are extended, and the larger the 
value of 0®, the more the disk emission will be suppressed, whereas the fraction of envelope 
obscured by the disk will be negligible. If the disk is seen close to edge-on, the envelope 
emission may eventually dominate. Figure 22 shows two of the most prominent objects, 
sources 177-341 (HST 1) and 182-413 (HST 10), as observed by with MAX on UKIRT 
and, for comparison, with the HST. These systems, with the disks seen edge-on, are clearly 
resolved at 10 /im with some evidence of a central peak. Size and orientation in the mid-IR 
correspond to that observed in the recombination lines. These images provide conclusive 
evidence that photo-evaporated dust can be heated in the photoionized globule and emit 
significantly in the mid-IR, with brightness comparable to that of the disk. It is clear that 
all possible combinations of angles 6a and 6§ may generate a variety of SEDs similar to 
those presented in Figure 16. All are characterized by an SED peaking at ~ 30 — 60 fim and 
relatively low fluxes at wavelengths shorter than 10 jum. 

Let us consider the impact of other basic parameters on the SEDs. In Figure 18, we 
present the SED obtained by placing our standard proplyd at various distances between 30" 
(d=0.065pc) and 150" (d=0.33pc) from the ionizing star. The SEDs are renormalized to the 
total stellar luminosity to compare their structure independently of the geometrical dilution 
of the signal with the distance. Still, the (relative) stellar flux in the visible increases with the 
distance, due to the corresponding reduction of optical depth within the globule described 
in Section 3.2 (see also Figure 1). In the far-IR the flux drops with the distance because the 
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outer disk is colder, as it receives less flux from the ionizing star, and it is therefore flatter, 
receiving less flux also from the central star. In comparison, the flux from the superheated 
layer is less affected by the distance from the ioninzing star, and for this reason the silicate 
peak in emission becomes relatively more prominent. 

Figure 19 shows the variation of the SED with Rout, the distance of the ionization front. 
We have considered R out — 100 AU, 200 AU and 300 AU. The emerging stellar flux is reduced 
due to the increase of optical depth in the envelope with the globule size (see again Figure 
1). The relative importance of the envelope emission increases significantly, but the SED 
is still dominated by the disk emission due to the special alignment of our standard case 
mentioned above. To facilitate the comparison with our standard case, we have neglected 
here the contribution from the hemisphere opposite to the ionizing star. 

The change of SED with the inner disk radius, Oj, is shown in Figure 20. The gap in 
the inner disk produces the expected decrease of the SED between 2 and 10 jum. Finally, 
the variation of SEDs with the stellar temperature is illustrated in Figure 21. Since we have 
kept the stellar radius constant, this plot also illustrates the variation of SEDs with the 
luminosity of the central star. Reducing the stellar temperature the disk emission becomes 
more and more important because the external contribution has been assumed to remain 
constant. 

A direct application of our model is presented in Figure 17, showing a fit to the available 
photometric data of source 177-341=HST1 (see Figure 22). We have assumed 9 d = 45° and 
Oq) = 75°, to account for a disk seen nearly edge-on. We used a distance d = 43'.'2, twice 
the projected distance from ^ 1 Ori C. Hillenbrand (1997) gives a spectral type later than 
K6. We have adopted a Ml spectral type with T* = 3750 K and radius R* = 2.18 R©. 
Concerning the geometry, we have assumed an ionization front at R out = 300 AU and a disk 
with an inner radius a, = 0.3 AU. To match the stellar photometry, we assumed A v = 2.2. 
The model predicts tr—1A5, so that approximately 0.6 magnitudes of visual extinction must 
be attributed to foreground extinction. In this fit, the envelope and disk have comparable 
brightness, the disk emission been provided by the disk interior. 

Given the number of parameters left unconstrained by the available observational data, 
it is not surprising that with our model we are able to fit a single photometric point at 10 /im. 
On the other hand, the 10 /im photometry cannot be readily explained with conventional 
disk models. Our model also provides a strong prediction for high mid- and far-IR fluxes 
from these systems. Future high spatial resolution observations at these wavelengths will 
asses the validity of our assumptions and provide insights into the evolution of circumstellar 
disks in the most typical environmental conditions. 
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A number of assumptions deserve some additional comment. First, we have assumed 
that the disk atmosphere remains well defined and treatable by the CG97 approximation even 
when the flaring angles become extreme (up to ~ 45°). In presence of photoevaporation, 
however, the conditions of hydrostatic equilibrium used to set the disk profile are no longer 
valid, and the locus of the disk "atmosphere" should be recalculated taking into account 
the radially diverging flow. Second, although we have a basic understanding of the globule 
physics, it is still very hard to constrain the actual structure from the observations, deriving 
information such as the density distribution. Third, the grain properties are critical, not 
only in what concerns the absorption efficiencies (Section 3.4), but also the albedo, that we 
neglected. Scattering processes affect the near- and mid-IR emission of disks seen nearly 
edge on. Also, it has been shown that dust scattering in an optically thin envelope can 
contribute to the increase of disk temperature at large distances from the central star (Natta 
1993), therefore adding a further contribution to the energy input at the origin of the disk 
flaring. Finally, we assumed in Section 2 a simplified model for the Orion Nebula. It is 
known (O'Dell 2001a, O'Dell 2001b) that the density in the Orion Nebula decreases with 
the distance from the interface with the Orion Molecular Cloud, where the main ionization 
front is located. Also, 9 1 Ori C is approximately three times closer to the main ionization 
front than to the foreground veil of neutral material, where a secondary ionization front is 
also located. The major effect of a more refined description would be a reduction of our 
estimates for the Lj-a energy density. We have seen that this is a secondary source of dust 
heating. A modulation over the r~ 2 dependency of the central radiation field would also be 
introduced. This effect may be significant in future studies aimed to the detailed comparison 
of the SED for different objects. 



7. Summary and conclusion 

We have explored the IR emission of circumstellar disks in the environment where star 
formation most typically occurs, i.e. a HII region powered by massive OB stars. This scenario 
applies in particular to the Orion Nebula, where the interaction of circumstellar disk with 
the environment has been directly resolved by the HST. We have build our model considering 
the four types of radiation relevant to the dust heating in a HII region, e.g ionizing (EUV) 
and non-ionizing (FUV) flux from the exciting star (f^Ori C), resonant Lj-a radiation and 
the remaining nebular radiation. We evaluated how these radiation sources affect the IR 
SED arising from: 

1. A spherical, homogeneous, optically thin, circumstellar globule, photoevaporated by 
the UV radiation of the HII region exciting star. The globule has a neutral core of 
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uniform density, and the photoevaporated atmosphere is treated following the Dyson 
(1968) model. With our assumptions, the most important parameter constraining the 
IR emission, the radial optical depth to the UV radiation, depends only on the distance 
from the ionizing star and on the globule size. The thermal emission peaks in the range 
10 — 30 /im, depending mostly on the dust composition. 

2. A non-accreting disk directly exposed to the nebular radiation. The disk is in hy- 
drostatic and radiative equilibrium, and treated following the prescription of Chiang 
and Goldreich (1997). We modify the CG97 scheme to account for the fact that the 
nebular radiation hits the disk surface with large angles. We follow the propagation 
of the various radiative fluxes through the disk superheated atmosphere, deriving the 
vertical temperature profile of the disk atmosphere. The disk faces receive unequal 
amount of radiation and present different flaring angles. In particular, the flux emitted 
from the face opposite to the ionizing star provides a model to the IR emission of the 
dark silhouette disks observed in the Orion Nebula. 

3. A combined system composed by the disk and the photoevaporated envelope, i.e. a 
photoionized proplyd of the type observed in the immediate surroundings of ^ 1 Ori C 
in the Orion Nebula. 

Depending on the distance and on the tilt angle of the disk with respect to the ionizing star, 
the disk flaring may be substantially higher than in the case of isolated disks. The tilt angle 
with respect to the Earth plays also a major role by hiding the central parts of the disk. The 
relative intensity of the disk vs. envelope emission varies with the tilt angle to the direction 
of the Earth. The high temperatures reached by the dust either at the disk atmosphere or 
within the envelope produce a SED peaking at 30 — 60 /im. 

We explore the dependency of the SEDs upon the tilt angle with respect to the Earth, 
the distance from the ionizing star, the size on the envelope, the inner disk radius, and the 
temperature of the disk's central star. The resulting SEDs are characterized by a broad 
peak of emission at 30-60 fim and are in general significantly different from those of isolated 
disks in low-mass star forming regions like Taurus-Auriga. Our model indicates that in the 
presence of an external radiation field, relatively evolved "Class 2" objects may display a 
SED peaking at mid-IR and far-IR wavelengths. The model explains the strong mid-IR 
excess we have recently detected on several sources in a 10 survey of the Orion Nebula. 
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Fig. 1. — Optical depth r° = a d riiRi as a function of the distance d from 9 1 Oii C. Six 
different globules are considered, with ionization front at R — 50, 100, 200, 300, 400, and 
500 AU. 
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Fig. 2. — Sketch of our assumed globule geometry with the radiative fluxes relevant for the 
dust heating. 
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Fig. 3. — Propagation of the flux from ^Ori C within the ionized wind. The two solid 
lines represent the FUV flux (upper line) and the EUV flux (lower line). The dashed line 
represents the dust optical depth, whereas the dot-dashed line represents the gas optical 
depth. The dotted line represents the ionization fraction. 
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Fig. 4. — Globule geometry. The left side shows our naming conventions for the optical 
depths, the right side for the distances. There is a simple relation between optical depth and 
distance (Equation 18), and corresponding quantities share the same index. Note, however, 
that t corresponds to the distance B4 — r. 




Fig. 5. — Attenuation of an isotropic radiation field within a homogeneous spherical globule 
with radial optical depth tr. Fom top to the bottom: 7^=0.01 , O.f , f , fO, 30. 
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Fig. 6. — Flux density along the globule axis. The vertical dotted line indicates the position 
of the ionization front at Ri = 100 AU. The fluxes represented are: non-ionizing FUV 
radiation from 6 |1 Ori C (solid thin line), ionizing EUV radiation from ^Ori C (solid thick 
line); radiation from the globule central star (long-dashed line), recombination line radiation 
from the nebular environment and from the ionization front (dot-dashed line); Ly-a radiation 
from the nebular environment and from the ionization front (short-dashed line), this last close 
to zero and mostly visible with a peak at the ionization front. 
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Fig. 7. — Temperature distribution along the globule axis resulting from the radiation sources 
represented in Figure 6. 
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log X(yLim) 

Fig. 8. — Spectrum emitted by a dusty globule with our standard set of parameters. The 
two solid curves are relative to type-1 grains (left) and type-2 grains (right). Each curve is 
the sum of the contributions from the neutral (dashed line) and ionized (dotted line) parts 
of the envelope. Self absorpion has not been considered. The stellar spectrum is not plotted. 
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Fig. 9. — Radius a d at which the disk heating from the central star equals that from the 
HII region as a function of the disk tilt angle Od and of the distance from ^Ori C. From 
top to bottom, the lines represent decreasing distances from 6 |1 Ori C. The dashed part of 
each line represents the case of pure nebular heating, occurring on the rear face of the disk 
(dark silhouette) within the HII region. The dotted lines at 9 d = 80° — 90° are relative to the 
case of a disk aligned almost edge-on with respect to the direction of ^Ori C. The stellar 
parameters are T* = 4, 500 K and = 2.5L . 
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Fig. 10. — Vertical sections of the disk. The continuous line represents the locus of the 
disk photosphere when both the flux from ^Ori C and the nebular radiation are taken into 
account. The dot-dashed line has been obtained by neglecting the flux from ^Ori C and 
is therefore representative of the dark face of the disk. The dotted line has been obtained 
taking into account only the radiation from the disk central star (as in CG97). 
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Fig. 11. — Temperature of the superheated layer in absence of external luminosity at dis- 
tances (top to bottom) a=3, 10, 30, 100 AU from the central star. 
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Fig. 12. — SED for a hydrostatic, radiative equilibrium disk directly exposed to the flux 
from ^Ori C. Solid line: stellar photosphere; dotted line: superheated layer; dashed line: 
inner disk. The thick solid line represents the total SED. The standard set of parameters 
with Draine and Lee (1984) silicates has been assumed. 
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Fig. 13. — Same as Figure 12, for a hydrostatic, radiative equilibrium disk directly exposed 
to the nebular radiation only. 
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Fig. 14. — SED of a flaring disk embedded within a photoevaporated envelope. Solid line: 
stellar photosphere; dotted line: superheated layer; dashed line: inner disk; dot-dashed line: 
outer envelope; thick solid line: total SED. 



-44- 




Fig. 15. — Image at 10 /im of the ^Ori C field. The arrows mark the position of ^Ori C. 
Source SC3 is the bright knot at located 1" to the West of 9 1 Ori C. The field size is 
approximately 15" x 15" and the pixel size is 0.26"/pixel. The image has been obtained with 
the MAX camera on the UKIRT telescope and is part of a large mosaic of the Orion nebula 
core (Robberto et al 2002, in preparation). 




Fig. 16. — SED of a flaring disk embedded within a photoevaporated envelope for different 
values of the tilt angle with respect to the Earth 9 e . The thick line refers to the standard 
case. 
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Fig. 17. — Photometry of the proplyd source 177-341 = HST 1 (crosses) and model SED 
assuming a flaring disk embedded within a photoevaporated envelope. See text for the fit 
parameters. The 10 /zm photometric point derives from the image presented in Figure 22. 
PHotometry at shorter wavelengths is from Hillnebrand (1997) and from Robberto et al. 
(2002). 
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Fig. 18. — SED of a flaring disk embedded within a photoevaporated envelope for different 
values of the distance from ^Ori C. The thick line refers to the standard case. 
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Fig. 19. — SED of a flaring disk embedded within a photoevaporated envelope for different 
values of the ionization radius R ou t- 
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Fig. 20. — SED of a flaring disk embedded within a photoevaporated envelope for different 
values of the disk inner radius a,-. 
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Fig. 21. — SED of a flaring disk embedded within a photoevaporated envelope for different 
values of the effective stellar temperature T*. 
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Fig. 22.— The prominent proplyd sources 177-341 (HST 1, left) and 182-413 (HST 10, 
right) imaged by the HST (top) and MAX on UKIRT (bottom). The HST images refer to a 
combination of recombination lines, namely [Nil] (red), Ha (green) and [OIII] (blue). The 
MAX images have been taken in N-broad band with integration time of the order of 600 
seconds on source (see Robberto et al. 2002 for details) 



